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O ' Abstract 

>>: 

Quark lepton universality inherent in grand unified theories based on SO (10) gauge group generi- 
cally leads to hierarchical neutrino masses. We propose a specific ansatz for the structure of Yukawa 
matrices in SO( 10) models which differ from this generic expectations and lead to quasi degener- 
ate neutrinos through the type-I seesaw mechanism. Consistency of this ansatz is demonstrated 
through a detailed fits to fermion masses and mixing angles all of which can be explained with 



a: . 

reasonable accuracy in a model which uses the Higgs fields transforming as 10, 120 and 126 rep- 



resentations of 50 (10). The proposed ansatz is shown to follow from an extended model based 
on the three generations of the vector like fermions and an 0(3) flavour symmetry. Successful 
numerical fits are also discussed in earlier proposed models which used combination of the type-I 
and type-II seesaw mechanisms for obtaining quasi degenerate neutrinos. Large neutrino mixing 
angles emerge as a consequence of neutrino mass degeneracy in both these cases. 
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I. INTRODUCTION 

Experiments over the years have revealed that 

(1) Two of the neutrino mixing angles are large as opposed to the small quark mixing angles. 

(2) Neutrino mass hierarchy is milder compared to quarks and extreme case of all neutrinos 
being quasi degenerate is still an allowed possibility 

Several independent reasons have been advanced [lH3| to understand feature (1) of the 
fermion spectrum but it may be that its answer lies in (2). Large mixing angles become quite 
natural if neutrinos are almost degenerate. They remain undefined in the exact degenerate 
limit. A small perturbation which leads to differences in neutrino masses can also stabilize 
all or some of the mixing angles to large values. Thus theory which predicts quasi degen- 
eracy has built in mechanism to explain large mixing angles. We present an SO (10) based 
unified description of fermion masses and mixing leading to hierarchical charged fermions 
and quasi degenerate neutrino masses. Large neutrino mixing emerges as a consequence of 
neutrino mass degeneracy. 

50(10) models provide a natural framework for understanding neutrino masses because 
of the seesaw mechanisms [l] inherent in them. Neutrino masses arise in these models from 
two separate sources either from the vacuum expectation value of the left-handed triplet 
(type-II) or from the right handed triplet (type-I) Higgs. It was pointed out [4, |5j long ago 
that combination of these two sources provides an interesting framework for understanding 
quasi degeneracy of neutrinos. In this approach, some flavour symmetry leads to degenerate 
type-II contribution and its breaking in the Dirac neutrino masses then leads to departure 
from degeneracy through the type-I contribution. This is realizable if type-II contribution 
dominates over the type-I which is not always the case[6j, |7j. Alternative possibility is that 
both degeneracy and its breaking arise from a single source namely type-I seesaw mechanism. 
This however requires a peculiar structure for the RH neutrino mass matrix Mr. It has 
been pointed out that the required structure can arise from the "Dirac screening" [8| or 
more generally from the application of the Minimal Flavour Violation [9| hypothesis to the 



leptonic sector jlO 



While these possibilities are known there does not exist a detailed study of all fermion 
masses and mixing in the context of realistic SO (10) models with quasi degenerate neutrinos 
and we address this question using (A) type-I mechanism alone and (B) combination of type-I 
and type-II mechanisms. 

We use supersymmetric SO (10) as our basic framework. Fermion masses arise in renor- 
malizable SO (10) models through their couplings to Higgs fields transforming as 10,126 and 
120 representations. One needs at least two of these fields to get fermion mixing and the 
minimal model with 10 and 126 has attracted a lot of attention|2j, la, LD, lu|. There have been 



studies of models with additional 120 also 12-141. In our context, we find that all the three 



Higgs representations are needed to obtain satisfactory fits to fermion masses and mixing. 
Starting with a supersymmetric 5*0(10), an effective Minimal Supersymmetric Standard 
Model (MSSM) is obtained by assuming fine tuning which keeps only two Higgs doublets 



lights. The final fermion mass matrices obtained after SO (10) and SU(2)l x U(l) breaking 
can be parametrized as |13|, Il4j : 

M d = H + F + G , 
M u = r (H + sF + t u G), 
M h = H-3F+ t t G } 
M D = r(H-3s F + t D G), 
M L = r L F, 
Mr = r-^ F, (1) 

where the matrices H and F are complex symmetric and G is an anti-symmetric matrix 
in generation space. H, F, G arise from the fermionic Yukawa couplings to the 10, 126, 120 
Higgs fields respectively, r, s, t u , £/, t D , r L , r R are complex parameters determined by the 
ratios of vacuum expectation values (vevs) and mixing among various Higgs doublets. The 
light neutrino mass matrix is given by 

M v = r L F - T R M D F~ x Ml = Mi 1 + Mi ■ (2) 

It is known that the above fermion mass structure allows different mixing patterns for 
quarks and neutrinos if type-II seesaw mechanism dominates [2|, lll| . Consider the limit in 
which the contribution of the 10-plet H dominates. In this limit, all the charged fermions 
are diagonalized by the same matrix and the Cabibbo-Kobayashi-Maskawa (CKM) matrix 
becomes proportional to identity. In the same limit, neutrino mixing with the type-II dom- 
inance is governed by F in eq.([2]) leading to non-trivial leptonic mixing. In fact, if only H 
dominates the charged fermion masses then one can obtain b-r unification which in turn 
drives the large atmospheric mixing & The existing fits QQ to fermion masses and mix- 
ing with type-II dominance are for the hierarchical neutrino masses. Degenerate neutrino 
spectrum can be obtained in this approach with an additional assumption: 

F = col (3) 

/ denoting a 3 x 3 identity matrix. The sub-dominant type-I contribution can then lead to 
the quark mixing and neutrino mass differences. 

The realization of the attractive type-II dominated scenario was found difficult in the 
context of the minimal model p, lZ|. It was found that parameter space favored by the 
overall fit to fermion masses suppresses the type-II contribution compared to the type-I. 
This motivates us to study degenerate neutrinos in the context of a purely type-I seesaw 
mechanism. A general framework based on flavour symmetry arguments to obtain quasi 



degenerate neutrinos in type-I seesaw was recently discussed [10] . It was found that the 
right handed neutrino mass matrix in an effective theory invariant under a spontaneously 
broken flavour symmetry Gf = 0(3); x 0(3) e x 0(3);, x U(1)r has a dominant term 

M R ^M T D M D + ..., (4) 



where each 0(3)/ acts on the flavour spaces of the lepton doublet (/ = /), the right-handed 
electrons (/ = e) and the right handed neutrino (/ = v) fields respectively and U(1)r 
corresponds to the right handed lepton number. The above form for Mr leads to the 
degenerate neutrinos through the type-I seesaw mechanism. An equivalent description at 
the SO (10) level can be obtained by imposing the following ansatz: 



aH 2 



(5) 



In the following, we will work out the detailed consequences of this ansatz for fermion masses 
and mixing and show that such an ansatz can be realized using a flavour symmetry. Since 
H is a symmetric matrix it can be diagonalized by a unitary matrix. 



U T H U = D 



H 



(6) 



where Dr- is a diagonal matrix with real elements. Without loss of generality, we can express 
the mass matrices in (jTJ) in an 5*0(10) basis with a diagonal H. This basis are obtained 
from eq. ([!]) by the replacement H — > Dh and 



U T FU = a(U T HU U^U* U T HU) = a D H V*D H . 



(7) 



G retains its antisymmetric form and we use the same notation for it and for various mass 
matrices in the new basis. From now on, we will work in this rotated basis. V = U T U in 
eq.Q is a symmetric unitary matrix which can be parametrized as 



V = PR^)U 12 (e,a)R 23 (4>)P 



with 



#23(0) 



(l 

cos' 
y sin < 
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- cos< 



/ 



U 12 (9,a) 



I cos 9 sin 6 \ 
sin 9 — cos 9 
e ia 



\ 



(8) 



(9) 



and P = Diag(e iai ,e ta2 , 1) is a diagonal phase matrix (One phase in P is absorbed in the 
complex parameter a in eq.([7|)). 

Before we present the detailed fits let us look at the implications of the ansatz eq.Q 
qualitatively. 

• Correct b-r unification and second generation masses are obtained if dominant con- 
tribution to the charged fermion masses come from the 10-plet, i.e. from H with a 
sub-dominant contribution from 126, 120 fields. Retaining only the H contribution, 
the ansatz, eq.([5]) implies that 



M[ = -r R M D F~ l M T D 



r 2 r R 



V + ..., 



(10) 



where the ... terms arise from the 126 and 120 contribution to the Dirac mass matrix 
Mjj- CKM matrix is unity in this limit while the neutrino mixing is determined from 



V. Diagonalization of V leads [15] to #23 = (ft, 9\i = § and #13 = where the angles 
9{j are angles defined in the standard parametrization of the leptonic mixing matrix in 
which #12 drives the solar and #23 the atmospheric neutrino oscillations. Thus ansatz 
in eq.([5]) can lead to correct description of the quark and leptonic mixing angles to 
zeroth order without requiring the type-II dominance as is commonly done. 

• If if in the original basis was real then V entering eq.(j7J) would be unity. In this 
case, all the fermion mixing vanish in the absence of the 120 contribution. Thus 
complex couplings and CP violation proves to be important in understanding large 
neutrino mixing within this approach. Numerically, we find that even after including 
120 contribution, one cannot get the correct mixing pattern with a real H. 

The mixing angles obtained at zeroth order with H dominance get corrected by the con- 
tributions from 126 and 120-plets. They induce non-zero quark mixing angles and perturb 
eq.flrOD: 

2 
Mi(M x ) = - T -^(V -GsaDH + toiGD^V -VD^G)) + 0{s\t 2 D )) , (11) 

The above neutrino mass matrix corresponds to an effective dimension five operator induced 
after integration of the right handed neutrino fields. Assuming that the heavy mass scale 
is close to the GUT scale and neglecting the effect of the Dirac neutrino couplings in the 
renormalization group (RG) evolution one can obtain [l| the low scale neutrino mass matrix 
as follows. Define the neutrino mass matrix M. v f(M x ) in the flavour basis as 

Mu f (M x ) = UjMu(M x )U? , (12) 

where TJ\ diagonalizes the charged lepton mass matrix M\ . The radiatively corrected neutrino 
mass matrix is then given by |l| 

Mvf(M z ) = I T M u {M x )ll , (13) 

where I T ~ Diag(l, 1,1 + e T ) and e r ~ "^liS?' 11 ^ 1, More detailed treatment then 
presented here would need to include RH threshold and ruining between the GUT and the 
RH mass scale etc. 



II. NUMERICAL FITS: TYPE-I SEESAW 

We now discuss detailed fits to fermion masses and mixing based on the ansatz ([5]) and 
the fermion mass matrices, eq.([T]). The latter are defined at the GUT scale M x , We use 
as our input the quark and lepton masses obtained at M x in the MSSM for tan/3 = 10, 



Msusy = ITeV and M GUT = 2 x 10 GeV [161 ] . The quark mixing angles do not appreciably 



change compared to the low scale values and we use the specific values used in earlier analysis 



presented in 13j,|l4|. Unlike in the previous works, the neutrino mixing angles are susceptible 



to change by the renormalization group (RG) evolution due to quasi degenerate nature of 
neutrinos. We include this effect as follows. Using the charged lepton mass matrix at Mx, 
we numerically determine the neutrino mass matrix in the flavour basis at Mx through 
eq. (fl~2l . Neglecting running of the mixing angles in Ui, the low scale neutrino mass matrix 
in eq.f lTHj) is numerically determined and used to obtain the observable neutrino masses and 
mixing angles. For neutrino masses and lepton mixings, we use the updated low energy 
values given in [17J. All our input data (central values and la errors) are summarized in 
Table dU). 



m d [MeV] 


1.03 ±0.41 


Am 2 oZ [eV 2 ] 


(7.59 ±0.20) x 10" 5 


m s [MeV] 


19.6 ±5.2 


Am 2 im [eV 2 ] 


(2.51 ±0.12) x 10~ 3 


m b [MeV] 


1063.6+^ 5 4 


sin 8f 2 


0.2243 ±0.0016 


m u [MeV] 


0.45 ±0.15 


sin0f 3 


0.0351 ±0.0013 


m c [MeV] 


210.3273+™ 


sin 9f 3 


0.0032 ± 0.0005 


m t [MeV] 


oo/lqq q+30267.6 

oz^oo.o_ 14768 g 


sin 2 9[ 2 


0.32 ±0.016 


m e [MeV] 


0.3585 ± 0.0003 


sin 2 # 23 


0.45 ± 0.092 


m M [MeV] 


75 671 ^+ - 0578 
'°- D ' io -0.0501 


sin 2 6[ 3 


< 0.049(3a) 


m r [MeV] 


1292.2±JJ 


$CKM 


60° ± 14° 



TABLE I: Input data used in this work (See text for details and references). 



We do the x 2 fitting to check the viability of the model as previously done in [7|, [Lj, Il4j . 
In this case we have total 25 real parameters (3 in Dh, 5 in V, 6 in G, real r, complex 
s, a, t u , ti, tr> ) which are fitted over 16 observables (9 charged fermion masses, 4 CKM pa- 
rameters, 2 leptonic mixing angles and Am 2 so J Am 2 atm ). Lepton mixings and Am^/Am 



2 
atm 



are independent of the overall neutrino mass (mo 



\TRT- 



■|) appearing in eq.( TTT|) , This thus 



remains an arbitrary parameter and the overall degenerate mass scale cannot be fixed. We 
also set r = — and minimize \ 2 with respect to the remaining 24 parameters. The results 
of the minimization are displayed in Table (ILTjl. Three of the viable solutions are displayed. 
We obtained the best fit value of x 2 — 2.038 for which all the observables are fitted within 
< 0.9cr. Solution(3) is also acceptable which fits all observables within < 0.7<r with excep- 
tion of down quark mass m^. The fits obtained here are better than the one obtained by 
Bertolini et al [7| in case of the minimal model with complete type-II seesaw dominance and 
hierarchical neutrinos. Unlike here, their fits have several observables which are > 1 or 2a 
away from the central values. 

We give below values of input parameters determined from \ 2 minimization in case of 
solution(l) of TabletfTll). 



D H = Diag. (0.00033658, 0.0149966, -0.757334) 



(14) 







Sol. 1 


Sol. 1 


Sol. 2 


Sol. 2 


Sol. 3 


Sol. 3 


No. 


Observables 


Fitted value 


Pull 


Fitted value 


Pull 


Fitted value 


Pull 


1 


m d [MeV] 


0.653677 


-0.917861 


0.678809 


-0.856564 


0.207819 


-2.00532 


2 


m s [MeV] 


17.5885 


-0.386821 


22.8346 


0.622041 


21.6923 


0.402361 


3 


m b [GeV] 


1.11131 


0.418721 


0.94463 


-1.0440 


1.05832 


-0.046348 


4 


m u [MeV] 


0.462718 


0.0847896 


0.461582 


0.0772103 


0.450825 


0.00549932 


5 


m c [GeV] 


0.210603 


0.0136849 


0.212603 


0.113153 


0.211727 


0.0695654 


6 


m 4 [GeV] 


63.6891 


-0.832404 


56.7159 


-1.14208 


67.6155 


-0.658038 


7 


ra e [MeV] 


0.358503 


0.00969691 


0.358516 


0.0525082 


0.358506 


0.0206782 


8 


m M [MeV] 


75.6719 


0.00734514 


75.6765 


0.0923818 


75.6711 


-0.0083064 


9 


m r [GeV] 


1.29219 


-0.00814429 


1.2921 


-0.0804718 


1.29223 


0.0218404 


10 


Am ltm 


0.0303514 


0.050109 


0.0302197 


-0.00862837 


0.0303237 


0.0377877 


11 


sin 9\ 2 


0.224205 


-0.0592102 


0.224193 


-0.0666865 


0.224306 


0.00359473 


12 


sin ^ 3 


0.0351308 


0.023704 


0.0347491 


-0.269936 


0.0350426 


-0.0441173 


13 


sin 0f 3 


0.00319336 


-0.0132867 


0.00322056 


0.0411106 


0.00315871 


-0.0825897 


14 


sin 2 6[ 2 


0.319801 


-0.0619079 


0.319568 


-0.076109 


0.321124 


0.0187774 


15 


sin 2 ^ 3 


0.481942 


0.313909 


0.437263 


-0.169787 


0.436492 


-0.178126 


16 


sin 2 8[ 3 


0.0195266 


- 


0.0463404 


- 


0.00288176 


- 


17 


5c KM [°] 


67.7227 


0.247333 


49.8678 


-0.422669 


56.4935 


-0.134071 


18 


5pA/ATs[°] 


53.9786 


- 


-52.7788 


- 


-66.9939 


- 




x 2 




2.038 




3.844 




4.684 



TABLE II: Three best fit solutions for fermion masses and mixing obtained assuming the type-I 
seesaw dominance. Various observables and their pulls obtained at the minimum are shown (See 
text for details). The bold faced quantities are predictions of the respective solutions. 

G 12 = (0.00167 + 0.00154J) ; G 13 = (0.0108 + 0.01017) ; G 23 = (0.191 + 0.0337) . (15) 



a = -(5.24 - 9.34/)GeV" 1 ; t t = 0.15 - 1.057; t u = 0.551 - 0.0847; 
s = (-3.21 + 3.327) x 10~ 4 - 



50.65°; # = -56.9°; a 



,t D = -(2.56 + 0.437) x TO -4 ; 
0.20°; ai = -39.55°; a 2 = 83.79° 



(16) 



Elements of Dh and G are expressed in GeV units. Note that <fi and 0/2 respectively 
determine the atmospheric and the solar mixing angles at Mx in the absence of perturbation. 
These values get stabilized at Mx to ~ 40.4° ; ~ —61.0° once the perturbations from 
126, 120 couplings are added. RG running changes them to the required values displayed in 
Table flLTj). 613 has not been included in our definition of \ 2 and its initial value was zero. 
This becomes non-zero but remains small in all the three solutions displayed. However, 
almost the entire allowed range in #13 is compatible with reasonable fits to other fermion 
masses as shown by the three solutions. Unlike the solar and atmospheric mixing angles, the 



ratio AmyAm^ im changes appreciably from 0.09 at M x to 0.03 at Mz by the RG effects. 
All these solutions predict large CP violating leptonic phase. The heaviest RH neutrino 
mass scale following from the above numerical fit would be approximately (see eq.((T]and 
ansatz (jSJ)) 

M 3 ttr^lalml « — ml ^2 x 10 13 GeV ( — — 
m V m o 

2 

where m = ^y-- Thus the RH neutrino mass falls below the GUT scale for this particular 
solution. 

Let us now illustrate how the ansatz (jHJ) can be obtained in a model from a flavour 
symmetry. A simple flavour symmetry to be used is 0(3) under which three generation of 
the 16-plet ip transform as triplets. The 0(3) breaking is introduced through a complex 
flavon field r\ transforming as spin 2. We need to introduce three generations of vector-like 
multiplets ^y + ^ v transforming as (16,3) + (16,3) under 50(10) x 0(3) and a U(l) x 
symmetry in order to realize eq.([SJ). The X-charges of [ip, \I/y,\I/y, rj, </>io, (fesr) are criosen 
respectively as (x, y, —y, l/2(y — x), —(x + y), —2y) with x ^ y. The general super potential 
invariant under 5O(10) x 0(3) x U(l)x can be written as: 

A A' 

W = Mm v m v + ^ v V V (J) m + 7*y^0io + — ^W^ + TT Tr tf^ + ( 17 ) 

Nip Nip 

The 0(3) and U(l)x breaking originates in the above super potential only from the Planck 
scale effects through the vev of the flavon field rj. The last two terms are the only terms 
which determine both the 10 and 126 Yukawa couplings once the heavy vector like fields are 
integrated out. The dotted terms correspond to terms suppressed by Mp. Here, the mass 
M of the vector like pair and the scale of the vev of rj lie above the GUT scale. The effective 
theory after integration of the vector like field is represented by 

w "< x wk"^" 1 "^ + mW p ^° • < 18 > 

where 

5' 

£ab = vlb + jTrr] 2 5 ab 

and a,b = 1,2,3 refer to the 0(3) index. This effective super-potential is also 50(10) x 
0(3) x U(l)x invariant. The Yukawa coupling H is proportional to the (£) and is a general 
complex symmetric matrix. The F is related to the square of H and satisfies the ansatz 
in eg. (151). The coupling to the 120 field can be generated by introducing a flavon field x 
with the U(l)x charge-2 x and transforming as a triplet of 0(3). This leads to the Yukawa 
coupling matrix G through the coupling 

^w P n ™ 

The above example thus illustrates how ansatz such as eq.([5]) can be realized. A detailed 
model along this line will require study of the details of the vacuum structure of the potential 
involving T], x and possibly additional fields for generating the right structure of the Yukawa 
couplings H, G. 

8 



III. NUMERICAL FITS: TYPE-II SEESAW 

We now turn to the numerical discussion of the ansatz (ED in which the contribution 



of 126 to fermion masses is assumed to be 0(3) invariant. The 0(3) breaking arises from 
the H and G contributions which lead to departure from degeneracy through the type-I 



seesaw. We shall not specify how this breaking occurs 18]. Such an ansatz for the type-II 
contribution was considered [5] in the specific context of SO (10). Detailed fits to fermion 
masses with recent data are however not presented in these works. We do this essentially 
following the same procedure as the one adopted in purely type-I case. We assume H, G to 
have the most general form. One could choose to work in a basis with a diagonal H. In this 
basis, eq.(j3]) gets changed to 

F = c V , (19) 

where V is a unitary symmetric matrix defined in eq.(|5J). In this basis, the charged fermion 
mass matrices can be obtained from eq.([T]) by replacing H with diagonal Dh, and F with 
CqV. The neutrino mass matrix, eq.([2]) can be written in the same basis as 

M v = mo (V - e M D V*Ml) (20) 

The parameter e controls the contribution from type-I seesaw which induces splittings in 
neutrino masses. 

We use these equations to fit all the fermion masses and mixing using the previous 
procedure. Results corresponding to the minimal case are displayed in Table fillip . The best 
fit solution we obtained here corresponds to x 2 = 6.0 which is acceptable for 16 data points 
from statistical point of view and all the observables except TO(, and m s are fitted with less 
than la accuracy. The obtained fit in the type-II case is however not as good as in the case 
of pure type-I seesaw combined with the ansatz (j5j). 

Numerical fits lead to e ~ 2 x 10~ 6 GeV~ 2 . Since the scale of Mp is set by the top mass the 
type-I contribution relative to the type II is given by em 2 ~ 10 -2 and type II contribution 
dominates as assumed. Now the overall scale of the RH neutrino mass is given by (see eq. (P 
and ansatz (J3j)) 



1 „ -„i K „ „ /0.2eV 



M, « -^— w 2 x 10 15 GeV 



is close to the GUT scale unlike the minimal models with type-II dominance but hierarchical 
neutrinos [y, |7[ . 

IV. SUMMARY 

Obtaining a unified description of vastly different patterns of quark and lepton spectrum is 
a challenging task. This becomes more so if neutrinos are quasi degenerate. We have shown 
here that it is indeed possible to obtain such a description starting from the fermionic mass 
structure, eq.((T]) that can arise in a general 5O(10) model. We considered two distinct 



No. 


Observables 


Fitted value 


Pull 


1 


m d [MeV] 


0.868041 


-0.395023 


2 


m s [MeV] 


12.2829 


-1.40714 


3 


m 6 [GeV] 


1.25634 


1.69141 


4 


m„[MeV] 


0.450489 


0.0032611 


5 


m c [GeV] 


0.210393 


0.00324503 


6 


m t [GeV] 


102.325 


0.883371 


7 


m e [MeV] 


0.358502 


0.00503107 


8 


m^MeV] 


75.6709 


-0.0111809 


9 


m T [GeV] 


1.29217 


-0.0244576 


10 


Am L 

Am atm 


0.0302538 


0.00659421 


11 


sin 0^ 2 


0.224154 


-0.0913125 


12 


sin 0f 3 


0.0351436 


0.033571 


13 


sin0? 3 


0.00326199 


0.123983 


14 


sin 2 0^ 2 


0.321168 


0.0214673 


15 


sin 2 23 


0.439779 


-0.14255 


16 


sin 2 6>[ 3 


0.0356836 


- 


17 


fcKM[°] 


49.7146 


-0.429864 


18 


<$PMiV,s[ ] 


-25.3338 


- 




x 2 




6.0 



TABLE III: The best fit solution for fermion masses and mixing obtained assuming the type-II 
seesaw dominance. Various observables and their pulls obtained at the minimum are shown (See 
text for details). 

possibilities based on purely type-I and the other based on the mixture of type-I and type-II 
seesaw mechanisms. Both these possibilities can lead to quasi degenerate spectrum if they 
are supplemented respectively with ansatz (j5]) and ()3]). We have shown through the detailed 
numerical analysis that these ansatz are capable of explaining the entire fermionic spectrum 
and not just the quasi degenerate neutrinos. Moreover, the origin of large leptonic mixing 
here is linked to the quasi degenerate structure determined by the matrix V providing yet 
another reason why quark and leptonic mixing angles are so different in spite of underlying 
unified mass structure. 
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